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STRESS  INTENSITY  FACTORS  IN  THE  THIRD-STAGE 
FAN  DISK  OF  THE  TURBINE  ENGINE 


INTRODUCTION 

Several  failures  of  the  TF-30  turbojet  engine  used  in  the  F-14  aircraft  have  been  attri- 
buted to  the  initiation  and  growth  of  fatigue  cracks  in  the  thi.d-stage  fan  disk  at  the  dovetail 
region  shown  in  Fig  1.  It  is  generally  believed  that  these  cracks  are  associated  with  an  inter- 
mittent high-frequency  engine  resonance  and/or  load  fluctuations  due  to  engine  power  or  speed 
changes.  Since  the  operational  static  load  ts  already  high,  only  moderate  fluctuations  are 
required  for  crack  initiation  followed  by  fatigue  growth.  In  either  case,  a reduction  in  the  static 
stress  level  extends  fatigue  life.  To  this  end,  an  experimental  stress  analysis  of  the  dovetail  was 
performed  by  Parks  and  Sanford'"  using  two-dimensional  photoelastic  models.  The  photoelastic 
analysis  showed  a stress  concentration  factor  of  5.2  in  the  disk  fillet  The  high  concentration  of 
tensile  stress  in  the  fillet  plus  the  fretting  action  at  the  edge  of  the  bearing  surface  are  the  two 
primary  factors  affecting  crack  initiation. 


Fig.  1 — TF-30  third-stage  blade/disk 
dovetail  region 
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Fatigue  cracking  occurs  in  a reproducible  manner  in  numerous  lugs  on  each  disk.  A typi- 
cal disk  removed  from  service  was  carefully  examined  at  NRL,  and  cracks  were  observed  in  15 
of  the  36  lugs.  In  each  instance,  the  crack  initiates  at  or  near  the  junction  of  the  fillet  and  bear- 
ing area  of  the  disk  lugs  as  illustrated  in  Fig.  2.  After  initiation,  the  crack  extends  upward 
across  the  lug  until  it  approaches  the  upper  fillet,  where  a compressive  stress  field  exists.  At 
this  point  the  crack  turns  through  an  angle  of  45  degrees  and  propagates  across  the  upper  region 
of  the  lug  as  shown  in  Fig.  3. 


\ 


Fig.  2 — Typical  fatigue  crack  in  a disk 
removed  from  service 


Fig.  3 — Final  phase  of  crack  extension 
after  the  crack  turns 


This  report  provides  a fracture  analysis  of  a crack  as  it  extends  from  the  initiation  point 
into  the  lug.  Of  particular  interest  is  the  stress-intensity  factor  or  the  strain  energy  release  rate 
at  the  crack  tip  as  the  crack  propagates  into  the  complex  stress  field  associated  with  the  dovetail 
region.  Knowing  the  stress  intensity  factor  along  the  path  of  the  crack  will  enable  a prediction 
to  be  made  of  the  propagation  life  of  the  fatigue  crack  from  the  Paris(2)  cyclic-growth  relation- 
ship 

da/dN  - C A Kn.  (1) 

where  da/dN  is  the  incremental  crack  extension  per  cycle,  A K is  the  range  in  stress  intensity 
factor,  and  n and  C are  material  constants. 

To  determine  the  stress  intensity  factor,  two-dimensional  photoelastic  models  3.409  times 
the  size  of  the  dovetail  region  of  the  disk  were  used.  Artificial  cracks  were  introduced  in  the 
models  by  sawing.  The  cracks  were  extended  incrementally  along  the  known  fracture  path,  and 
isochromatic  fringe  patterns  were  recorded  for  each  crack  length. 

The  isochromatic  patterns  in  the  neighborhood  of  the  crack  tip  indicated  that  the  stress 
field  was  of  the  mixed-mode  type  and  that  a far-field  stress  <rox  parallel  to  the  crack  was  of 
significant  magnitude.  As  a method  of  determining  the  opening  stress  intensity  factor  K,  and 
the  shearing-mode  stress  intensity  factor  K„  from  isochromatic  fringe  loops  under  these  condi- 
tions was  not  available,  it  was  necessary  to  develop  a data-analysis  procedure.  This  was  accom- 
plished by  using  a modified  near-field  solution  for  the  central  crack  problem  to  obtain  a non- 
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linear  equation  relating  K,,  K„,  and  <rox  to  the  fringe  order  N and  the  fringe  position  (r,  6) 
relative  to  the  crack  tip.  Solution  of  this  nonlinear  equation  and  the  extraction  of  Kh  K,h  and 
o-0x  was  obtained  using  the  Newton-Raphson  iterative  method  combined  with  a least-squares 
method  to  minimize  error  due  to  inaccuracies  in  measuring  fringe  parameters.  The  results 
obtained  for  K / and  Kn  were  accurate  to  10  to  15  percent  based  on  individual  determinations 
and  to  about  5 to  7 percent  based  on  mean-value  estimates  when  several  determinations  were 
made. 

Results  from  the  model  analysis  for  K,  and  K„  were  scaled  to  the  prototype  based  on  a 
centrifugal  load  of  24,900  lb  (111  kN)  applied  to  the  prototype.  An  effective  stress  intensity 
factor  Keff,  related  to  Kt  and  K/h  was  determined  as  a function  of  crack  length. 


EXPERIMENTAL  PROCEDURE 

The  experimental  approach  was  similar  to  the  approach  used  by  Parks  and  Sanford"’  in 
their  photoelastic  stress  analysis.  A complete  description  of  the  model  geometry,  the  loading 
fixture,  and  model  manufacturing  methods  is  contained  in  Ref.  1.  The  models  used  in  this 
study  were  made  from  the  same  templates  as  employed  by  Parks  and  Sanford  and  had  in-plane 
dimensions  3.409  times  the  prototype  dimensions. 

The  models  were  fabricated  from  Homalite  100*  sheets  0.255  in.  (6.48  mm)  thick.  A 
crack  was  simulated  in  the  model  by  using  a fine-pitched  jewelers  saw  0.018  in.  (0.5  mm)  thick 
with  40  teeth  per  inch.  The  crack  tip  was  sharpened  with  a swiss  jewelers  file’*  with  a radius  of 
about  0.008  in.  (0.2  mm). 


The  model  was  loaded  in  a sodium  diffused-light  polariscope,  and  isochromatic  fringe  pat- 
terns were  photographically  recorded.  The  crack  was  extended  incrementally  by  sawing  and 
filing;  then  the  loading  and  photographing  sequence  was  repeated.  A total  of  23  combinations 
of  load  and  crack  length  were  photographed. 


Typical  isochromatic  fringe  patterns  obtained  for  three  lengths  of  cracks  are  illustrated  in 
Fig.  4.  An  enlargement  of  the  fringe  loops  at  the  crack  tip  is  presented  in  Fig.  5.  Inspection  of 


Fig.  4 — Typical  isochromatic  fringes 
for  crack  extension  in  the  disk 
dovetail 


*G  & L Industries,  Wilmington,  Delaware. 
••Grobet  Vallorbe  No.  6 jewelers  file. 
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these  fringe  patterns  shows  that  the  loop  tilt  angle  6m  (defined  in  Fig.  6)  and  the  maximum 
radius  rm  associated  with  the  fringe  loops  above  and  below  the  crack  line  are  not  the  same. 
Thus  the  cracks  are  propagating  in  a mixed-mode  field  with  both  Kt  and  Kn  affecting  fatigue 
crack  extension. 

METHODS  OF  ANALYSIS 

A method  to  determine  mixed-mode  stress  intensity  factors  from  isochromatic  fringe 
loops  has  previously  been  developed  by  Smith  and  Smith'3*.  Their  method,  which  neglects  the 


Fig.  6 — Characteristic  shape  of  upper 
and  lower  loops 
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nonsingular  far-field  stress  <rox,  can  be  applied  only  when  the  axes  of  the  upper  and  lower 
fringe  loops  lie  along  a common  line  Examination  of  Fig.  5 indicates  that  this  condition  is  not 
satisfied  for  the  problem  being  considered  here.  For  this  reason  it  was  necessary  to  extend  the 
method  of  analysis  to  include  the  <rox  term. 


For  a crack  subjected  to  both  tensile  and  shearing  loads  the  stress  field  in  the  neighbor- 
hood of  the  track  tip  can  be  approximated <4)  as 
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where  the  coordinates  are  measured  from  the  crack  tip. 

The  maximum  in-plane  shear  stress  rm  is  related  to  the  Cartesian  components  of  stress  by 

(2rm)2  = (<ry  - crx)2  + (2t xy)  2.  (3) 


The  stress-optic  law  in  photoelasticity  which  relates  the  fringe  order  N to  the  maximum 
in-plane  shear  stress  is 

2 rm  = Nfjh,  (4) 

where  fa  is  the  material  fringe  value  and  h is  the  model  thickness. 

Substituting  Eqs.  (2)  and  (4)  into  (3)  gives 

( Nfjh )2  = — — [(AT/  sin  8 + K„  cos  0)2  + (K„  sin  0)2] 

Li t r 

+ p-JJ  sin  4 [Kt  sin  0(1  -I-  2 cos  8) 

•Jlirr  1 

+ A"//(l  + 2 cos2  0 -I-  cos  0)J  + (TgX.  (5) 


The  size  and  shape  of  the  isochromatic  fringe  loops  are  related  to  Kh  K„ , and  <rox 
through  Eq.  (5).  It  is  possible  to  use  this  relation  in  two  ways  — first,  to  assume  that  Kh  K„, 
and  <rox  are  known  and  to  use  Eq.  (5)  to  determine  the  geometry  of  fringe  loops  of  a given 
order  or,  second,  to  assume 'the  geometry  of  the  fringe  loop  is  known  and  to  use  Eq.  (5)  to 
determine  Kh  Kn , and  <rox. 

Thus  Eq.  (5)  forms  the  basis  of  all  of  the  methods  used  to  determine  Kh  Kn,  and  <rox 
from  isochromatic  fringe  patterns.  In  this  analysis  three  approaches  were  followed  in  develop- 
ing solutions.  In  each  approach  the  set  of  equations  to  be  solved  was  nonlinear  in  the  un- 
knowns Kf,  A',/,  and  <rox.  These  equations  were  sob  ed  using  the  Newton-Raphson  method  as 
applied  to  systems  of  equations15'.  Each  of  the  approaches  will  be  described. 


Classical  Approach 

In  the  first  determination  of  K,  and  aox  from  symmetrical  fringe  loops,  performed  by 
Irwin(6),  it  was  noted  that,  for  closed  fringe  loops  near  the  crack  tip,  a point  exists  where 
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This  point  was  assigned  coordinates  (rm,  9m).  This  relation  gave  Irwin  a second  equation  to  use 
with  a simplified  form  of  Eq.  (5)  (AT//  terms  eliminated)  and  enabled  him  to  solve  for  K,  and 
(t  ox  in  terms  of  the  fringe  order  Nand  the  position  parameters  rm  and  9 m. 

This  approach  can  also  be  used  in  the  mixed-mode  mode  analysis  by  using  information 
from  two  loops  in  the  near-field  of  the  crack.  Differentiating  Eq.  (5)  with  respect  to  9 , setting 
9 *=  0m  and  r = rm , and  using  Eq.  (6)  gives 

f(K„  K,„  <rox)  = ~ ~ (*,2  sin20m  + 4 K,  K„  cos2 9m  - 3 K,2,  sin2 9m) 

L-nrm 


2 (j.  Q 

— . °x  sin  ~r-  [Kl(cos9m  + 2cos20m)  - K„(2s\n29m  + sin0)] 
2 


, 1 

+ — cos  — - 
2 2 


[AT/(sin0m  + sin20m)  + K„(2  + cos20m  + cos0m)] 


= 0, 


(7) 


where  the  coordinates  rm  and  9m  define  the  point  on  any  fringe  loop  which  is  farthest  removed 
from  the  crack  tip,  as  illustrated  in  Fig.  6.  Rewriting  Eq.  (5)  in  a form  suitable  for  use  with  the 
Newton-Raphson  method  gives 


. / TS 


+ Z.-L.  sin  [K,  sin  9m  (1  + 2 cos  9m) 

y/2nrm  2 

+ K„(  1 + 2 COS2  9„  + COS  Ol  + <*lx  ~ (Nm  /<r/ft)2  = 0-  (8) 


where  Nm  is  the  fringe  order  corresponding  to  the  point  (rm,  9m). 

When  the  classical  approach  is  used,  the  stress  intensity  factors  K,  and  K„  and  the  remote 
stress  (Tox  are  determined  from  isochromatic  data  taken  from  two  points  — one  on  each  of  two 
fringes.  In  concept  any  two  loops  can  be  used;  however,  in  practice  more  accurate  results  are 
obtained  if  one  fringe  is  taken  from  the  set  of  loops  above  the  crack  line  and  the  other  fringe  is 
taken  from  the  set  of  loops  below  the  crack  line. 

Substituting  the  radii  rm  and  the  angles  9m  from  these  two  loops  into  Eq.  (7)  gives  two 
independent  relations  in  terms  of  the  parameters  of  interest.  The  third  equation  is  obtained  by 
using  Eq.  (8)  with  data  from  either  the  upper  or  the  lower  loop.  The  three  equations  obtained 
in  this  manner  are  of  the  form 


/Mi.  *//.  O = o. 

/,(*„  K„,  <tox)  = 0.  (9) 

gu(Ki.  K/i,  o’  ox ) = 0, 

where  the  subscripts  u and  / refer  to  upper  and  lower  loops  respectively. 

Although  Eqs.  (9)  can  be  solved  in  closed  form,  the  algebra  becomes  quite  involved,  and 
a simpler  approach  using  a numerical  procedure  based  on  the  Newton-Raphson  method  was 
employed.  To  describe  the  Newton-Raphson  method,  consider  an  arbitrary  function  hk  of  the 
form 

M*,.  *//.  «<J  - 0.  W 
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where  k = 1,2,  and  3.  If  initial  estimates  are  made  for  Kh  K,h  and  <tox  and  substituted  into 
Eq.  (10),  the  equality  is  not  satisfied,  since  the  estimates  are  in  error.  To  correct  the  estimates, 
a series  of  iterative  equations  based  on  a Taylor  series  expansion  of  hk  are  written  as 


(9*),+ 1 


(/!*),  + 


dhk 

dK, 


A K,  + 


BK„ 


MC„  + 


9o  ox 


(11) 


where  the  subscript  / refers  to  the  ith  iteration  step.  A K,,  AA„,  and  \vox  are  corrections  to 
the  previous  estimates.  The  corrections  are  determined  so  that  (hk)l+\  = 0,  and  thus  Eq.  (11) 
gives 

A = ~(hk),.  (12) 


dhk 


dK, 


aa,  + 


dhk 


9A„ 


A K„  + 


dhk 


dcr. 


Applying  Eqs.  (12)  to  Eqs.  (9)  and  solving  for  the  correction  terms  &K,,  A K,h  and  A cr0 
yields,  in  matrix  notation, 
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dfu 
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are  given  by 

<+i  = 

(K,),+ 

AA,, 

,+i  = 

(K„),+ 

i+i  = 

(o-0x),  + b<Tox 

-l 


(13) 


(14) 


The  convergence  of  this  method  is  rapid,  and  usually  four  or  five  iterations  are  sufficient 
for  obtaining  precise  estimates  of  K,,  K„,  and  <t ox.  A computer  program  written  in  BASIC, 
with  a file  name  MXDMD4,  was  prepared  to  compute  Kh  K„,  and  crox  according  to  Eqs.  (13) 
and  (14).  This  program  is  listed  in  Appendix  A. 


Selected-Line  Approach 

It  is  possible  to  simplify  Eq.  (5)  by  considering  only  the  isochromatic  data  along  the  line 
defined  by  9 — ir.  Then  Eq.  (5)  reduces  to 


WJh)2 


4 Kh  . AK„<r„  , , 

+ — + <T„ 


2 nr  y/2irr 
which  is  independent  of  K,.  Equation  (15)  can  be  written  as 

2 K„ 


Nfjh  - ± 


U-nr 
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(15) 


(16) 


7 


SANFORD  AND  DALLY 


The  negative  sign  is  selected,  because  Ktl  < 0 for  isochromatic  intercepts  with  the  6 = tt  line.* 
Consider  two  isochromatics  Nx  > jV2,  which  intercept  the  9 = tt  line  at  positions  rx  and  r2  as 
shown  in  Fig.  7 and  select  the  minus-sign  solution  in  Eq.  (16).  Then 


(V,  /„  2 K„ 

(T° * h 72777 

and 


(17a) 


Nif„ 


2 K, 


h 72777- 

Subtracting  Eq.  (17b)  from  Eq.  (17a)  to  eliminate  ctox  gives 

„ fa  77  V'l  r 


2 — 7^2 


(jV,  - tfj). 


(17b) 


(18) 


which  is  negative,  since  Nx  > N2  and  r2  > rx.  Finally  aox  is  obtained  from  Eq.  (iS)  and  Eq. 
(17a)  as 


To  determine  Kh  reevaluate  Eq.  (5)  with  9 = tt/2  and  obtain 

(Nfjh)2  = (AT,2  + AT,2,)  + ~=  (K,  + K„)  + o-2,. 

Equation  (20)  is  quadratic  in  terms  of  AT,;  thus 

—b  + vQ)2  — 4 ac 


K,= 


2 a 


(20) 


(21) 


*A  similar  development  can  be  made  when  Ac//  >0,  except  that  the  isochromatics  intersect  the  tower  edge 
of  the  crack  at  0 — — tt  . 
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- 


where 


a — 1/2-rr  r3, 

b - CToJyJlT  r3,  (22) 

c - (Kn/lirr  J + ( AT,/  cr  0x!  + <*lx  ~ (N}fjh)2. 

Note  that  N}  and  ri  are  defined  in  Fig  7.  The  plus  sign  is  selected  in  front  of  the  radical  in  Eq. 
(21),  since  K,  > 0. 

Least-Squares  Approach 

The  least-squares  approach  fits  Eq.  (5)  to  the  data  at  a number  of  arbitrary  points  over  the 
fringe  field  in  the  neighborhood  of  the  crack  tip  using  the  unknowns  Kh  K„,  and  a-nx  as  the 
fitting  parameters.  The  fitting  process  is  similar  to  that  used  in  the  classical  approach  in  that  the 
Newton-Raphson  method  is  employed;  however  the  least-squares  approach  uses  data  from  k 
arbitrary  points,  where  k > 3. 

Applying  the  iteration  scheme  suggested  in  Eq.  (12)  to  k equations  of  the  form  given  in 
Eq.  (8)  results  in  an  overdetermined  set  of  linear  equations  in  terms  of  the  unknowns  AA/( 
AAW,  and  Ao-0JI  of  the  form 


dft 

dft 

dft 

bK, 

bK„ 

do- ax 

ft 

dft 

dg2 

dft 

ft 

- - 

dAT, 

bK„ 

do  ax 

ft) 

\ 

dft 

dft 

dft 

bK, 

bK„ 

do  ax 

For  brevity  rewrite  the  matrices  in  Eq.  (23)  as 

U)  - [a]  IAAT1. 


AA, 

A*// 


(23) 


(24) 


It  can  be  shown  that  [A AT]  in  Eq.  (24)  can  be  determined  in  the  least-squares  sense  by  multiply- 
ing both  sides  from  the  left  by  la] r,  where  la]  ris  the  transpose  of  [a]: 

[a]rld  = la] r [a]  [A AT]  (25) 


or 

where 

and 

Finally, 

where  [c]'1  is  the  inverse  of  lc] . 


[ d ] = lc]  [A*], 

l d)  - lain*] 


lc]  = la] T la]. 
IAAT]  - Id"'  Id, 


(26) 


(27) 
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The  solution  of  Eq.  (27)  gives  AX/,  A K/h  and  Ac t ox,  which  are  used  to  correct  initial  esti- 
mates of  K,,  AT /,,  and  a ox  to  obtain  a better  fit  of  the  function  g to  the  k data  points.  By  use 
of  the  matrix  operations  available  in  BASIC,  a computer  program  with  a file  name  MXDMD5 
was  written  to  determine  K/,  K //(  and  c r ox  using  data  from  ten  points  over  the  fringe  field. 
This  program  is  listed  in  Appendix  B. 

VERIFICATION  OF  METHODS 

All  three  methods  were  verified  and  were  found  to  give  accurate  results  (to  ±0.1  per- 
cent), providing  exact  input  data  were  used.  The  data  required  for  verification  were  generated 
by  substituting  selected  values  of  K,,  Kn,  and  cr ox  into  Eq.  (5).  The  resulting  expression  was 
then  used  to  determine  the  fringe  order  A at  specified  ( r , 9)  pairs. 

Analysis  with  real  data  is  less  accurate,  because  experimental  errors  will  occur  in  deter- 
mining the  r and  9 associated  with  the  location  of  a point  on  an  isochromatic  fringe.  Errors  in  9 
of  ±2  degrees  and  in  /-of  ±0.005  in.  (0.13  mm)  should  be  expected.  The  errors  are  primarily 
due  to  difficulties  in  locating  the  crack  tip  at  the  end  of  the  saw  cut. 

The  classical  approach  is  perhaps  the  most  error  prone  of  the  three  methods,  since  9m 
must  be  determined  for  both  the  upper  and  lower  loops.  Considering  the  difficulty  in  locating 
the  crack  tip  and  in  locating  the  point  corresponding  to  rm  on  both  isochromatic  fringe  loops, 
errors  in  9 m of  ±5  degrees  may  be  anticipated.  An  example  was  run  with  computer  program 
MXDMD4  such  that  9m  on  the  upper  loop  was  increased  by  5 degrees  and  9m  on  the  lower 
loop  was  decreased  by  5 degrees.  The  error  in  the  predictions  of  K,,  Kn,  and  <tox  was  13.9, 
15.6,  and  38  percent  respectively. 

The  selected-line,  classical,  and  least-squares  methods  were  compared  by  using  the  analyt- 
ically generated  point  plot  of  the  fringe  pattern  shown  in  Fig.  8.  Measurements  of  r,  9,  and  N 
were  taken  from  this  point  plot  using  the  same  techniques  which  are  used  with  actual  photo- 
elastic fringe  patterns.  The  data  with  typical  measurement  errors  were  processed,  and  the 
results  of  the  error  analysis  are  shown  in  Table  1.  These  results  indicate  that  the  least-squares 
method  with  the  program  MXDMD5  provided  the  most  accurate  estimates  of  K , and  Kt,. 

ANALYSIS  OF  DATA 

Negatives  of  the  fringe  patterns  were  projected  in  a photographic  enlarger  onto  the  data 
sheet  shown  in  Fig.  9.  The  projected  size  of  the  image  was  approximately  10  times  the  actual 
model  size.  Five  data  points  associated  with  each  fringe  loop  were  marked  on  the  data  sheet. 
These  data  points  were  usually  at  the  extremity  of  a fringe  loop,  where  changes  in  r with  9 are 
small  and  measurement  error  is  minimized.  The  scale  distance  (5  mm)  from  the  negative  was 
marked  on  the  data  sheet. 

Measurements  were  made  with  proportional  dividers  to  determine  the  radius  associated 
with  each  data  point.  Values  of  angle  9,  radius  r,  and  fringe  order  N associated  with  ten  to  40 
data  points  were  made  from  each  negative  of  the  isochromatic  pattern. 

These  data  were  used  as  input  for  the  MXDMD-5  program  to  determine  Kh  Klh  and  <tox. 
Ten  data  points  were  used  for  each  determination,  five  from  an  upper  loop  and  five  from  a 
lower  loop.  When  more  than  ten  data  points  were  available,  different  combinations  of  data  sets 
were  employed  and  several  estimates  of  Kh  Kn,  and  <rox  were  made.  In  these  instances,  aver- 
ages of  the  individual  determinations  were  taken.  The  mean  values  of  Kh  Klt,  and  <rox  are 
presented  in  Table  2. 
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Table  1 — Comparison  of  the  three  Methods 
on  the  Basis  of  Typical  Measurement  Errors 


True 

Value 

Least-Squares  Method 

Selected-Line  Method 

Classical  Method 

Predicted 

Value 

Error  (%) 

Predicted 

value 

Error  (%) 

Predicted 

Value 

Error  (%) 

K(  (psi-/m.) 

318.974 

317.5 

0.5 

304.8 

4.5 

298.4 

6.4 

K|(  (psiVin.) 

-293.184 

-297.1 

1.3 

—308.2 

5.1 

-308.0 

5.1 

,r«v(Psi) 

354.865 

345.3 

2.7 

380.7 

7.2 

357.9 

0.8 
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Table  2 — Summary  of  Results  for  K , K u , and 
< t for  the  Photoelastic  Model 

OX 


Record  No. 

Number 

of 

Readings 

P 

Load 

(lb) 

Crack 

Length 

(in.) 

*/ 

(psi^/m.) 



K„ 

(psiVin) 

"ax 

(psi) 

Before  Crack  Turns 

D-2-6  & 7 

1 

800 

0.096 

573 

-569 

2482 

D-2-8  & 9 

1 

1000 

0.096 

693 

-706 

2974 

D-l-3 

1 

200 

0.122 

439 

-216 

327 

D-l-4 

1 

200 

0.185 

475 

-159 

238 

D-l-5  & 6 

5 

200 

0.241 

348 

-219 

377 

D-l-7  & 8 

6 

200 

0.357 

318 

-285 

357 

D-1-9&  10 

7 

200 

0.466 

512 

-203 

298 

D-l-11  & 12 

7 

200 

0.622 

539 

-281 

318 

D-l-13  & 14 

5 

200 

0.710 

549 

-221 

134 

D-l-15  & 16 

3 

200 

0.843 

392 

-334 

152 

D-l-17  & 18 

5 

200 

0.955 

351 

-444 

213 

After  Crack  Turns 

D-l-19  & 20 

1 

200 

0.100 

793 

-365 

370 

D-l-21  & 22 

3 

ISO 

0.190 

782 

-184 

431 

D-l-23  & 24 

1 

200 

0.190 

1038 

-198 

533 

D-l-25  & 26 

2 

100 

0.296 

815 

-86 

131 

D-l-27  & 28 

1 

100 

0.296 

786 

-163 

154 

D-l-29  & 30 

1 

100 

0.410 

1044 

244 

-326 

D- 1-31  & 32 

1 

100 

0.518 

1196 

347 

-519 

SCALING  FROM  MODEL  TO  PROTOTYPE 

Conversion  of  the  results  from  the  photoelastic  model  of  the  turbine  disk  to  results  which 
correspond  to  the  actual  prototype  requires  the  use  of  a scaling  relationship.  The  scaling  equa- 
tion for  K,  or  K,i  is 

( Kl.ll)p  = >A7  ^2  ^3 (K/,//)m-  (28) 

where  the  subscripts  p and  m refer  to  the  prototype  and  model  respectively.  The  terms  A,,  A2, 
and  A 3 in  Eq.  (28)  are  scaling  factors  defined  as  follows: 

A | = the  in-plane  scaling  factor  - wm/wp  - 3.409,  (29) 

where  wm  is  the  center-to-center  dimension  between  two  adjacent  blades  on  the  model  and  is 
3.808  in.  (96.7  mm)  and  wp  is  the  corresponding  dimension  on  the  prototype  and  is  1.117  in. 
(28.38  mm); 

A 2 - the  thickness  scaling  factor  - hm/hp  -0.111,  (30) 

where  hp  is  the  length  of  the  dovetail  slot  and  is  2.30  in.  (58.34  mm)  and  hm  is  the  model 
thickness  and  is  0.255  in.  (6.48  mm);  and 

A 3 = the  load  scaling  factor  - Pp/Pm  — 124.5,  (31) 

where  Pp  - 24,900  lb  (111  kN)  and  P„  - 200  lb  (890  N). 
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The  scaling  equation  for  the  length  a of  the  crack  is 

a„-am/\v  (32) 

These  scaling  relations  assume  that  the  crack  in  the  prototype  propagates  in  a two- 
dimensional  manner  (as  through-thickness  cracks  with  a straight  front).  Inspection  of  the  frac- 
ture surface  of  the  lugs  which  have  been  pulled  from  the  disk,  as  shown  in  Fig.  10,  indicate 
that  the  cracks  propagate  with  three-dimensional  characteristics.  The  cracks  initiate  at  a point 
and  do  not  extend  along  the  entire  thickness  of  the  disk.  Instead  the  crack  front  exhibits 
significant  curvature.  Because  of  these  differences  in  the  crack  shape  between  the  model  and 
the  prototype,  the  results  presented  here  should  be  considered  as  a first  approximation.  More 
exact  predictions  would  require  three-dimensional  photoelastic  studies  to  more  adequately 
model  the  crack  shape. 


Fig.  10  — Typical  fracture  surface  showing  the  three- 
dimensional  characteristics  of  a fatigue  crack 


Numerical  values  of  K,  and  K „ related  to  a load  on  the  prototype  of  Pp  - 24,900  lb  (111 
kN)  are  given  in  Table  3 as  a function  of  crack  length  ap.  Inspection  of  Table  3 shows  that 
both  K,  and  Ku  vary  between  narrow  limits  as  the  crack  grows  from  0 to  0.280  in.  (7.1  mm) 
deep.  To  assess  the  influence  of  the  combined  effect  of  K / and  Kth  it  is  appropriate  to  consider 
the  total  strain  energy  release  rate  G,  which  is  given  by 

G — G,  + G„,  (33) 

where,  for  plane  strain, 

G,  - K}  (34a) 

and 

G„  - Kn-  (34b) 

in  which  E is  the  modulus  of  elasticity  and  v is  Poisson’s  ratio. 

Substituting  Eqs.  (34)  into  Eq.  (33)  gives 

G - K}„.  (35) 
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Table  3 — Stress  Intensity  Factors  X/  and  K 
in  the  Third-Stage  Fan-Blade/Disk  P = 24,900  lb 

p 

Centrifugal  Load 


*/ 

(ksi  Vm.) 

Ku 

(ksi  Vin.) 

a 

(in.) 

Keff 

(ksiin.) 

Before  Crack  Turns 

3.65 

-3.63 

0.28 

5.15 

3.54 

-3.60 

0.028 

5.05 

11.20 

-5.51 

0.035 

12.48 

12.12 

-4.06 

0.054 

12.78 

8.88 

-5.59 

0.071 

10.49 

8.11 

-7.27 

0.105 

10.89 

13.06 

-5.18 

0.137 

14.05 

13.75 

-7.17 

0.182 

15.51 

14.01 

-5.64 

0.208 

15.10 

10.00 

-8.52 

0.247 

13.14 

8.96 

-11.33 

0.280 

14.44 

After  Crack  Turns 

20.23 

-9.31 

.029 

22.27 

26.60 

-6.26 

.056 

27.33 

41.59 

-4.39 

.087 

41.82 

40.11 

-8.31 

.087 

40.96 

53.28 

12.45 

.120 

54.71 

61.03 

17.70 

.151 

63.54 

where 

*,//  - V*/2  + *//■  (36) 

The  effective  stress  intensity  factor  Kef/,  which  governs  the  rate  of  crack  propagation 
da/dN  according  to  Eq.  (1),  increases  from  zero  to  a value  of  about  14  ksiVirv  as  the  crack 
extends  to  a depth  of  about  0.100  in.,  as  shown  in  Fig  11.  Continued  crack  extension  from 
0.100  to  0.280  in.  occurs  at  nearly  constant  Krff. 

After  extending  a distance  of  0.280  in.,  the  crack  approaches  the  compressive  stress  field 
near  the  upper  fillet  of  the  lug.  The  value  of  A,  decreases  and  |X,/|  increases  until  the  crack 
turns  through  an  angle  of  about  45  degrees.  After  turning,  X,  becomes  dominant  and  Kel, 
increases  markedly  with  further  crack  extension.  It  is  believed  that  the  rate  of  crack  propaga- 
tion will  increase  rapidly  after  the  crack  turns. 


ERROR  ESTIMATE 

In  most  of  the  cases  examined  here,  several  fringe  loops  were  available  for  analysis  and 
30  or  40  data  points  could  be  obtained  from  regions  on  the  fringe  loops  where  r was  not  chang- 
ing rapidly  with  respect  to  9.  In  these  instances,  the  data  could  be  grouped  into  different  sets  of 
ten  data  points,  and  multiple  determinations  of  the  X field  could  be  made.  An  example  of  the 
scatter  in  the  determination  of  X,,  K,h  and  <rm  is  shown  in  Table  4 for  a typical  set  of  data 

. 
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The  mean  value  x,  standard  deviation  Sx,  and  coefficient  of  variation  Sjx  were  deter- 
mined and  are  also  listed  in  Table  4.  These  statistical  parameters  indicate  that  errors  in  indivi- 
dual determinations  of  Kt  of  about  ±10  percent  and  in  Kn  of  about  ±13  percent  can  be 
expected  about  68  percent  of  the  time.  If  the  error  is  due  to  random  causes,  then  the  averaging 
process  improves  the  accuracy.  The  standard  deviation  of  the  mean  S L is 

S-  - SJfn, 

where  n is  the  number  of  measurements. 

The  values  of  S-  and  S-/x  listed  in  Table  4 illustrate  the  improvement  in  predictions 
based  on  mean  values.  In  this  case,  errors  of  ±3.6,  ±5.0,  and  ±5.7  percent  in  Kh  K„ , and 
(Tm  can  be  anticipated  68  percent  of  the  time.  Higher  errors  would  be  anticipated  32  percent  of 
the  time. 
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Appendix  A 
PROGRAM  MXDMD4 


00001  REM  A PROGRAM  EMPLOYING  THE  CLASSICAL  METHOD  OF  ANALYSIS 
0000?  REM  FOR  DETERMINING  KI.K2  AND  SIGMA  SUP  OX 

00003  REM  THE  NUMERICAL  ANALYSIS  IS  PASED  ON  THE  NFWTON  PAPHSON  METHOD 

00004  REM 

00009  REM  CAUTION  THETA  FOR  LOWER  LOOP  IS  NEGATIVE. 

00010  PRINT  "INPUT  THETA, RADIUS  AND  FR1NGF  ORDFR  FROM  UPPER  LOOP") 
00015  INPUT  T 1 , R 1 , N I 

00020  PRINT  "INPUT  THFTA, RADIUS  AND  FRINGE  ORDER  FROM  LOWER  LOOP" ! 
00025  INPUT  T2.R?,N? 

00030  PRINT  "INPUT  MATERIAL  FRINGE  VALUE  AND  MODFL  THICKNESS"! 

00035  INPUT  P, H 

00040  PRINT  "INPUT  INITIAL  FST1MATF  OF  K|"l 
00045  INPUT  K 1 
00050  K2-0 
0C055  K3=0 

00060  PR  I NT"NUMPER  OF  ITERATIONS  M IS"! 

00065  INPUT  M 

00070  REM  COMPUTE  TRIG  FUNCTIONS 
00075  T1 =3. 1 4| 59*T1 / 180 
00080  T2  = 3. I 41 59*T2/ 180 
00065  S0  = SGRT<?*3. I 41 SS*P I ) 

00090  SI  =SIN(TI > 

00095  S?=C0SCT1) 

00100  S3=SlNCTl/?> 

00105  S4=C0SCTl/2) 

0C1I0  S 5=  SORT <2*3. I 41 59*  R2  > 

00115  S6=SIN(T2> 

00120  S7=C0S(T2> 

00125  S8  = S.IN<  T2/2) 

00130  S9=C0S(T2/2) 

00135  RFM  COMPUTE  COEFFICIENTS  OF  TAYLOR’S  SFPIES. 

00140  A1 =SI**2/S0**2 

00145  A2=4*S1*S2/S0**2 

00150  A3=(4*S?**?*S1**?>/S0**? 

00155  A4=P*S3*SI*( l*2**2)/S0 

00160  A5=?*S3*( 1+?*S2**2*S2)XS0 

00165  AR=(NI*P/H)**2 

00170  PI sS6**?/S5**2 

00175  P2=4*S6*S7/S5**2 

00180  R3=<  4*S7**2  + S6**?)/S5**2 

00185  P4=?*S8*S6*< 1*2*S7)/S5 

00190  E5=2*S8*( 1 +2*S7**2*S7)/S5 

00195  PP s<N?*P/H)**2 

00200  Cl =2»S1*S2/S0**2 

00205  CP=4*(S?**?-SI**?)/S0**2 

00210  C3=-6*S1*S?/S0**2 

0021 5 C4=<?/S0)*(S3*«  S?*?*S2**2-24SI**2)*t . 5* $4) ♦ ( S 1 ♦?♦ S I * S?>  > 

00220  C5*(?/S0>*< -S3*<  4*5. |*S?*SI ) ♦< • 5* S4 ) * ( ?♦ SP**2- S I **2* S2) ) 

00??5  Dl »?*Sf*S7/S5**2 
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00230  D2=A*<S7**2-S6**2>/S5**2 
00235  D3=-6*S6*S7/S5**2 

002 AO  DA  = (2/SS>*<  S3*  C S7*2*S7**2-2*S6**2> ♦( . 5* S9 > * < S6*2*  S6* S7 > > 
002  A5  D5=<2/S5)*<-S8*<A*S6*S7*S6>+<.5*P9>*(2*S7**2-S6**2*S7>) 
00250  PRINT 

00255  PRINT  " K SUF  I K SUB  2 SIOMA  SUB  0X" 

00260  PRINT  KI.K2/K3 

00265  DIM  F<3)»IJ(3»3>»V<3»3>»Y<3»3)»X(3> 

00270  REM  FFl-IN  ITERATIVE  S0LUTI0N 
00275  FOP  L= 1 T0  M 

00230  F < I >=At*Kl**2+A2*Kl*K2+A3*K2**2*Aa*Kt*K3*A5*K2*K3+K3**2-AS 

0028  5 F ( 2 > =R 1 *K I **2*R2*K 1 *K2»B3*K2**2*P4*K I *K3* F 5*K2*K 3+K3**2- FS 

00290  F<  3>  = DI*K l**2  + D2*K 1*K2  + D3*K2**2*Da*K1 *K3  + DS*K2*K3 

00295  Ut 1» n=2*A|*KI*A2*K2+AA*K3 

00300  U< 1»2>=A2*KI +2*A3*K2+A5*K3 

00305  U< 1 » 3)=Aa*K I +A5*K2*2*K3 

00310  U<  2» I>=2*B1*K1*P2*K2+B4*K3 

00315  U(2»2)=82*KI+2*P3*K2*B5*K3 

00320  U<2# 3>=BA*K I *B5*K2+2*K3 

00325  UI3< 1 >=2*D1*K1 *D2*K2+D4*K3 

00330  U< 3.2) = D2*K I +2*D3*K2*D5*K3 

00335  U<3j3)*Da*K1*D5*K2 

003 AO  MAT  V=INV(U) 

003A5  MAT  Y = V*F 
00350  MAT  X»{- I )*Y 
00355  K|=K1*X(  I ) 

00360  K2=K2*X< 2 ) 

00365  K3=K3+X(3) 

00370  PRINT  KIJK2IK3 
00375  NEXT  L 
00380  END 
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Appendix  B 
PROGRAM  MXDMD5 


00001  RFM  A PROGRAM  WHICH  USES  THE  LEAST  SCUARFS  METHOD  TOGETHER  WITH  THE 

00002  REM  NEWTON  PAPHSON  METHOD  T0  DETERMINE  K1.K2  AND  SIGMA  SUE  OX 

00003  RFM  FROM  AN  OVER  DETERMINED  SYSTEM  OF  EGUATI ONS* 

00004  RFM  THE  PROGRAM  IS  WRITTEN  TO  ACCEPT  10  DATA  POINTS. 

00005  REM 

00010  PRINT  "INPUT  INITIAL  ESTIMATE  OF  K 1 " 

0001 5 INPUT  K I 
00020  K2  = 0 
00025  K 3 = 0 

00030  PRINT  "INPUT  DATA  FOR  THETA  <T>»  RADIUS  tR>  AND  FRINGE  ORDER  fN>" 

00035  FOR  J=1  TO  10 

00040  INPUT  T< J).R< J).Nt J> 

00045  Tt J>=3. 1 41 S9*TtJ>/ 180 
00050  NFXT  J 

00055  PRINT  "MATERIAL  FRINGE  VALUE  AND  MODEL  THICKNESS"! 

00060  INPUT  P»H4 

00065  RFM  CONVFRT  UNKNOWNS  TO  NORMAL  I ? FD  FORM. 

00070  K 1 =K I *H4/P 

00085  PRINT  "SPECIFY  NUMPFR  OF  ITERATIONS"! 

00090  INPUT  0 

00095  DIM  GC  I0>*  At  10.  3>#  V.C3.  10).C<3.  3>»D<3>»U<  3*  3>»  Vt  3/  3).  X(  3) 

00100  PRINT  " K SUB  1 K SUB  2 SIGMA  SUB  OX" 

00105  FOR  1=1  TO  10 
00110  SOt I >=2*3. 1 41 59*Rf I) 

00115  SI ( I > = SINC  T(  I > > 

00120  S2t I >=C0S< Tt I)) 

00125  S3C I > = SINtTt  I >/2> 

00130  S4( 1 >=C0StTt  I )/2> 

00135  A| C I ) = (S1 C I > >**2/S0C I > 

00140  A2( I > = 4*S1 < I )*S2< I > / SO  t 1 > 

00145  A3< I > = ( 4*t S2< 1 > >**2+<S1 ( I > >**2>/S0< I > 

00150  A4C I >=2*S3t I >*SI t I >*t 1 +2*S2t I > > / SORT < SOt I > > 

00155  A5< I >=2*S3< I >*< l+2*t  S2< I > >**2+S2t I >>/SCRTt SO  t I ) > 

00160  A6t l >=tNt 1 > >**2 
00165  NEXT  I 

00170  RFM  BEGIN  ITERATIVE  SOLUTION. 

00175  FOP  L= 1 Te  G 

001R0  E I =0 

00185  FeR  1=1  TP  10 

00190  FtI>=A5tI>*K2*K3+K3**2-AftI> 

00195  G< 1 > = A » II )*K1**2+A2t 1 >*K 1*K?+A3t 1 >*K2**2*A4t  I)*K1*K3+E<1) 

00200  At I > | >=2*A1 ( I )*K 1 ♦ A2 1 I >*K2*A4t I )*K3 
00205  At  I * 2 > = A2 C I )*KI+?*A3< 1 ) *K2*A5t I > *K3 
00210  At I«3)=A4( I >*K1*A5C I >*K2*2*K3 
00215  E 1 *E 1 ♦ Gt I ) **2 
00220  NFXT  I 
00225  MAT  W=TPNtA> 

00230  MAT  C * W * A 
00235  MAT  0*  W*G 
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00240  MAT  U= I N V (C) 

00245  MAT  V=(- I >*U 
00250  MAT  X*  V*D 
00255  Kl =KI + X< 1 ) 

00260  K2  = K2*XC 2) 

00265  K3*K3*X(3) 

00270  PRINT  K1*P/H4IK2*P/H4JK34F/H4 
00275  NEXT  L 

00280  PRINT  "RELATIVE  ERROR- 
0028  5 PRINT  El 
00290  END 


